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Abstract 



The temporal and spatial coherence widths of the microscopic process by 
which a neutrino is detected are incorporated in the quantum mechanical 
wave packet treatment of neutrino oscillations, confirming the observation of 
Kiers, Nussinov and Weiss that an accurate measurement of the energies of 
the particles participating in the detection process can increase the coherence 
length. However, the wave packet treatment presented here shows that the 
coherence length has an upper bound, determined by the neutrino energy and 
the mass-squared difference, beyond which the coherence of the oscillation 
process is lost. 

PACS number: 14.60.Pq 
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A complete understanding of neutrino oscillations must take into account the localization 
of microscopic processes by which a neutrino is produced and detected. This localization is 
appropriately described by a wave packet treatment of neutrino oscillations [0-0 (different 
treatments are discussed in |10|]ll|| ). As the authors of |J noticed, in the quantum mechanical 
wave packet approach presented in j|] the dependence of the oscillation probability on the 
temporal and spatial coherence widths of the detection process was neglected. In this brief 
report, we wish to incorporate, in a simple and straightforward way, the temporal and spatial 
coherence widths of the detection process in the quantum mechanical wave packet description 
of neutrino oscillations and show that, as an immediate consequence of this, performing an 
accurate measurement of the energies of the particles involved in the detection process leads 
to an increase of the coherence length for neutrino oscillations, as was noticed for the first 
time in |J. 

Let us consider a neutrino of flavor a produced by a weak interaction process at the 
origin of the space-time coordinates and detected at a distance L after a time TQ by a 
weak interaction process capable of detecting a neutrino of flavor (3. As in [|J, we describe 
the neutrino propagating between the production and detection processes with the time- 
dependent state in the Schrodinger picture 
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where E a (p) = \Jp 2 + rr? a . This state is a superposition of mass eigenstate wave packet 
states (labeled by the index a) weighted by the complex-conjugated elements of the mixing 
matrix U of the neutrino fields. For simplicity, we consider only one space dimension in 
the source-detector direction and we assume that the mass eigenstate wave functions in 
momentum space ip a (p;p a , a p p) have the following gaussian form, which enables us to carry 
out an analytical calculation of the transition probability: 
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where p a are the average momenta of the different mass eigenstates and a p p is the momentum 
width of the wave packets. The average momenta p a are determined by the kinematics of the 
production process taking into account the masses m a of the mass eigenstates. We assume 
that all the mass eigenstates are extremely relativistic, i.e. p a 3> m a . In this case AE ~ Ap 
and a p p is determined by the minimum between the energy and momentum uncertainties 
of the production process (a possible dependence of the momentum widths from the index a 
can be estimated to be very small and negligible for relativistic neutrinos). Hence, from the 
uncertainty principle it is clear that a p p is determined by the maximum between the spatial 
and temporal coherence widths of the production process. 



1 Taking into account the spatial and temporal coherence widths of the production and detection 
processes, x = 0, t = and x = L, t = T are their average space-time coordinates. 
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We assume that the gaussian wave functions (0) are sharply peaked around the corre- 
sponding average momentum, i.e. a p p <C E^/m a , with E a = E a (p a ). Under this condition 
the energy E a (p) can be approximated by E a (p) ~ E a + v a (p — p a ), where v a = p a / E a is the 
group velocity of each wave packet. In this case, the neutrino wave function in coordinate 
space = (x\v a (t)) is easily calculated after a gaussian integration to be given by 

Va) , (3) 

with the orthonormal states \ v a ) belonging to the mass Fock space ({v a Wb) = 8 a b). The width 
<j x p of the mass eigenstate wave packets in coordinate space is related to the momentum 
width a v p by the uncertainty relation u x p<J p p = 1/2. From the considerations presented 
after Eq.(||D it is clear that the value of a x p is given by the maximum between the spatial 
and temporal coherence widths of the production process. 

Let us consider now the detection process which takes place at a distance L and after a 
time T from the origin of the space-time coordinates and is capable of detecting a neutrino 
with flavor (3. In Q the amplitude of the flavor changing process was obtained by pro- 
jecting the state |z/ Q (L,T)), which describes the propagating neutrino, on the flavor state 
\vp) = J^aUpaWa) ■ This procedure neglects the temporal and spatial coherence widths of 
the detection process. In order to take them into account, the detected neutrino must be 
described by a wave packet state analogous to that in Eq.(|l]): 

\ v Pi = ^2 U k dp^ a (p;p a ,a pD ) \u a (p)) . (4) 

a J 

The value of the momentum width a p p> is given, for relativistic neutrinos, by the minimum 
between the energy and momentum uncertainties of the detection process. The state \vp) 
does not evolve in time because it does not represent a propagating neutrino. The average 
momenta p a in the mass eigenstate gaussian wave functions describing the detected neu- 
trino are the same as those of the corresponding mass eigenstate gaussian wave functions 
describing the neutrino propagating between production and detection. This property is a 
consequence of causality: after the average momenta p a have been determined by the kine- 
matics of the production process, the mass eigenstates propagate between the two processes 
and determine the kinematics in the detection process. For example, the moduli of the av- 
erage momenta p a of the mass eigenstate wave packets of a muon neutrino produced in pion 
decay at rest are fixed by the kinematics of the process. When this neutrino is detected, for 
example, by scattering with a nucleus at rest, each mass eigenstate determines with its aver- 
age momentum p a and its mass m a a different value for the momenta of the recoil particles. 
Neutrino oscillations are observed if the different mass eigenstates are detected coherently, 
i.e. if the differences of the energies and momenta of the recoil particles corresponding to 
different mass eigenstates are smaller than the energy and momentum uncertainty of the 
detection process. 

Taking into account that the detection process takes place at a distance L from the origin 
of the coordinates, the wave function in coordinate space describing the detected neutrino 
is given by 
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where a x o is defined by the uncertainty relation a x o <J p d = 1/2. Hence, the value of a x o is 
given by the maximum between the spatial and temporal coherence widths of the detection 
process. 

The amplitude of the flavor changing process is given by the overlap 



A af3 (L,T) = / dx(vp(x - L)\v a (x, T)) . 



The integral over x is gaussian and leads to 
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with 



(8) 



This relation is very important in that it clearly shows that the width a x which determines 
the coherence of the flavor changing process depends on the spatial and temporal coherence 
widths of both the production and detection processes. The amplitude (0) has the same 
form as that in [|J, with an important difference: the width a x in [|J was determined only 
by the production process, whereas in Eq.(^) also the detection process is properly taken 
into account. 

In order to calculate the oscillation probability, the mass eigenstate energies can be 
approximated by 
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where E is the energy determined by the kinematics of the production process for a massless 
neutrino and £ is a dimensionless quantity that is determined by energy-momentum conser- 
vation in the production process to first order in m 2 a / E 2 . The quantity £ is typically of order 
unity; for example, in neutrino production by pion decay at rest we have £ ~ 0.2. With this 
approximation we have p a ~ E — (1 — £) m 2 a /2E and v a ~ 1 — m 2 /2E 2 . From Eqs.(^) and 
(0), the probability of v a — > va transitions is given by 
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In principle, P a p(L,T) is a measurable quantity, but in all realistic experiments the 
distance L is a fixed and known quantity, whereas the time T is not measured and can have 
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any value, because the source and the detector typically operate for times much longer than 
the oscillation times 4nE/ Am 2 b . Therefore, the quantity that is measured in all experiments 
is the oscillation probability P a p(L) at a fixed distance L, given by the time average of 
P a p(L, T). After integrating over T and imposing the normalization condition Y^r Papi.L) = 
1, we obtain 



P a p{L) = Ka U (3aU a bUp b exp 
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where 
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are the oscillation wavelengths and coherence lengths, respectively, and 
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Equation (|TT|) contains, in addition to the usual expression for the neutrino oscillation prob- 
ability (the first term in the exponential with F a b = 1), two additional factors, the first being 
the second term in the exponential that takes into account the coherence of the contributions 
of different mass eigenstate wave packets and the second being the additional factor F a t. The 
factor F a b is equal to unity if a x <C \L°^\, which is a necessary condition for the observation 
of neutrino oscillations that must be satisfied by any realistic experiment. If this condition 
is not satisfied, the interference among different mass eigenstate wave packets is washed out 
and only a constant transition probability P a p = J2 a \U aa \ 2 \Up a \ 2 can be observed. 

It is important to notice that the integration over the oscillation time is a crucial step in 
the wave packet treatment of neutrino oscillations in space, because it allows us to eliminate 
the time degree of freedom, which is not measured in realistic experiments. This elimination 
is done at the classical level by integrating the probability, which is a classical quantity, and 
without any unphysical assumption on the energies and momenta of the mass eigenstates, 
which are determined by the production process 
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An unphysical assumption is, for 
example, the imposition of equal energy to the different mass eigenstates, i.e. £ = 0, which 
would eliminate the time dependence of the oscillatory term in the probability (PH|). 
From Eqs.([D and (|T2|) one can see that the coherence length ' r " 1 ' 
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L c ° b is proportional to 
Hence, L^ h is dominated by the largest among the temporal and spatial 



coherence widths of the production and detection process. In particular, a precise determi- 
nation of the energies (or momenta) of the particles involved in the detection process implies 
a small a P £> and a large a x o, leading to a large coherence length L c a ° h h ~ A\f2a xD E 2 / |A 



,.. - , m lb\ 

(if a p o cTpp) . Hence, as noted in || , the coherence length can be increased by measuring 

accurately the energies of the particles participating in the detection process. However, even 
if, at least in principle, the coherence length can be increased without limit by an extremely 
precise measurement of the energies of the detection (or production) particles, the oscilla- 
tions lose the coherence for a x > L°^ c , because the factor F ao becomes important and has 
the effect of suppressing the interference of u a and v h for a x ^> L° a s b c . Physically this is due to 
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the fact that the spatial or temporal width of the detection (or production) process becomes 
larger than the oscillation length, leading to the washout of the oscillations|. In order to 
understand the reason of this washing-out let us consider, for example, the case in which 
the detection process has a small spatial coherence width and a large temporal coherence 
width that dominates the total coherence width a x . As shown by Eq. (|i0|) , in this case the 
interference between the mass eigenstates a and b is not suppressed if the average detection 
time T does not differ from L/v a and L/vb by more than ~ a x . This is due to the fact that 
the modes of the detection process corresponding to the different mass eigenstates oscillate 
coherently during a time interval of width a x around T, even when the mass eigenstate wave 
packets are far from ZQ. Then it is clear that if a x is larger that the oscillation length, the 
phase of the interference term depends crucially on the value of T. In this case, if T is not 
measured the average of the transition probability over T washes out the interference. 
Therefore, using Eq . fll2|) one can see that the coherence length has the upper bound 

i c 2 p3 

LT < E {LZf = (14) 

beyond which the coherence length loses its meaning, because the coherence of the process 
is lost for any value of the distances L. 

In conclusion, we have shown that the temporal and spatial coherence widths of the 
detection process can easily be incorporated in the quantum mechanical wave packet treat- 
ment of neutrino oscillations. As a result, we confirm the observation presented in |J that 
an accurate measurement of the energies (or momenta) of the particles participating in the 
detection process can increase the coherence length. However, the wave packet treatment 
presented here shows that the coherence length cannot be increased beyond the upper bound 
given by Eq.([L4]) without losing the coherence of the oscillation process. 



2 Notice that a similar effect is obtained if the distance L from the source is not known and the 
time-dependent probability P a p(T) is obtained by averaging P a p{L,T) over L. Such a situation is 
realized, for example, in the calculation of the effects of neutrino oscillations in the early universe. 

3 Notice that if o x d is dominated by the temporal coherence width of the detection process, 
the gaussian approximation in (||) for the wave function of the detected neutrino introduces an 
unphysical interaction for times smaller than the time of arrival of the propagating neutrino to 
the detector. In this case a more realistic approximation can be obtained by inserting in (||) and 
@ a 6{x — L + g x d), where a X D is the spatial coherence width of the detection process. This 
approximation requires a numerical solution of the integral in (|^) and will be discussed in detail 
elsewhere 
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